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ABSTRACT 



Aims. Study the contradictory magnetic field strength distributions retrieved from independent analyses of spectropolarimetric obser- 
vations in the near-infrared (1.56 yum) and in the visible (630 nm) at internetwork regions. 

Methods. In order to solve this apparent controversy, we present simultaneous and co-spatial 1.56 yum and 630 nm observations of 
an internetwork area. The properties of the circular and linear polarization signals, as well as the Stokes V area and amplitude asym- 
metries, are discussed. As a complement, inversion techniques are also used to infer the physical parameters of the solar atmosphere. 
As a first step, the infrared and visible observations are analysed separately to check their compatibility. Finally, the simultaneous 
inversion of the two data sets is performed. 

Results. The magnetic flux densities retrieved from the individual analysis of the infrared and visible data sets are strongly correlated. 
The polarity of the Stokes V profiles is the same at co-spatial pixels in both wavelength ranges. This indicates that both 1.56 yum 
and 630 nm observations trace the same magnetic structures on the solar surface. The simultaneous inversion of the two pairs of 
lines reveals an internetwork full of sub-kG structures that fill only 2% of the resolution element. A correlation is found between the 
magnetic field strength and the continuum intensity: equipartition fields (S ~ 500 G) tend to be located in dark intergranular lanes, 
whereas weaker field structures are found inside granules. The most probable unsigned magnetic flux density is lOMx/cm^. The 
net magnetic flux density in the whole field of view is nearly zero. This means that both polarities cancel out almost exactly in our 
observed internetwork area. 
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1. Introduction 

Solar magnetic fields as seen in magnetograms form a par- 
ticular pattern outside active regions, in the so-called quiet 
Sun. The largest polarimetric signals are confined to the 
supergranular boundaries forming the photospheric network. 
Inside these large-scale convective cells it is possible to 
detect magnetic fields with a more disperse character and 
with smaller polarimetric signals. These areas are called the 
internetwork. The magnetic structures of the network can 
be modeled either with a flux tube model or with the 
MISMA hypothesis (MIcroStructured Magnetic Atmosphere, 
ISanchez Almeida & Landi DegFInnocentil Il996h . Both models 
reproduce the observed Stokes profiles and agree on the exis- 
tence of nearly vertical magnetic fields of kG strength occupy- 
ing only som e 10-20 % of the reso l ution e lement. Ob servations 
at 1.56 um (iMudach & SolankiL [19921: iLiiiL fT99l. 5 25 nrn 
JStenflol [T973 b and 630 nm ("S anchez Almeida & Litei. l2000l) 
give consistent results. Unfortunately, this is not the case for the 
internetwork, with recent studies based on high quality spectro- 
polarimetric data at 1 .56 yum and 630 nm yielding contradictory 
results. 

The first observation of Stokes V spectra of an internet- 
work region was carried out by Keller et al. (1994) . These au- 
thors observed the pair of Fe i lines at 525 nm and analysed 



four average Stokes V sp ectra using the line ratio technique 
(iHoward & Stenflg, Il97lh . Their conclusion was that the mag- 
netic field strength was below 1 kG (500 G) with a probability 
of 95% (68%). Later, Lin ( 1995) measured directly the Zeeman 
splitting of the Fe i line pair at 1 .56 fim to infer the magnetic 
field strength. He showed that the typical magnetic field strength 
in the internetwork was well below kG. The retrieved kG struc- 
tures were correlated with network patches. A few years later, 
Lin & Rimmele (1999), using the same pair of infrared lines, 
deduced field strengths between 0.2 and 1 kG. 

The success of the Advanced Stokes Polarimeter dLites et al.L 
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1993) filled the literature with works concerning the internet- 
work magnetism using the Fe i pair of lines at 630 nm. Most 
of these studies obtained that the field strength distribution 
at the internetwork presents a peak at kG. These kG fields, 
even if they occupy a small fraction of the Sun's surface (0.1- 
1%), contained most of the magnetic flux of the internetwork. 
Applying a Principal Component Analy sis (PCA) inversion, 
ISocas-Navarro & Sanchez Almeidal (l2002h also recovered mag- 
netic field strengths of the order of kG. The database for this 
PCA algorithm consisted of synthetic profiles generated un- 
der the MISMA hypothesis. The most surprising fact is that, 
anal ysing the sam e data in terms of a Milne-Eddington inver- 
sion, iLitesI (I2OO2I) obtained fi eld strengths clearly below 1 kG . 
Using the line ratio technique. 'Domin guez Cerdena et al.](l2003h 
confirmed the presence of ubiquitous kG fields in the internet- 
work. They used speckle reconstructed data in order to achieve 
a spatial resolution of 0.5". 
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iKhomenko et alj (l2003h carried out the study of an inter- 
network region at 1 .56 yum, measuring the Zeeman splitting to 
deduce the magnetic field strength. These authors obtained that 
the internetwork magnetic field strength was well below the kG 
regime, with a probability distribution function that could be re- 
produced by a decreasing exponential with an e-folding factor of 
250 G. 

From all these works, one can summarize that studies us- 
ing the 630 nm or the 1.5 fim lines have led to contradictory 
results: the infrared lines indicate fields below 1 kG, while the 
visible lines indicate a dominance of kG fields. Two plausi- 
ble explanations have been presented to explain such a discrep- 
ancy. One is based on the idea that visible and near-infrared ob- 
servations contai n incompatible information about the internet- 
work magnetism dSocas-Navarrol 120031 : ISanchez Almeida et al.L 
l2003h . The other suggests that the line pair at 630 nm is not 
reliable for retrie ving the magnetic field strength at the internet- 
work (.Bellot Rubio & Colladosl 120031: iMartinez Gonzalez et al.L 
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l2006tlL6pez Ariste et al.Ll2006h . 

ISocas-Navarrol (|2003|) suggested that, if kG and sub-kG 
fields co-exist in the same resolution element, the visible lines 
would be sensitive to the kG contribution, whereas the infrared 
lines would trace the sub-kG structures. The explanation comes 
from the different sensitivi ty of the spectral lines t o the mag- 
netic field. According to Bellot Rubio & Colladosl 1*2003), this 
may be the situation in the case of having kG fields occupying ~ 
20% of the resolution element and sub-kG fields filling the rest. 
If the per centages change, the premise s eems not valid anymore. 
However. lSocas-Navarro & LitesI (|2004|) found observational ev- 
idence indicating the coexistence of sub-kG and kG fields in the 
internetwork. Their analysis was based on a three-component in- 
version, two of the components being magnetic with fixed field 
strengths of 1700 and 500 G and the other one being field-free. 
The observations were compatible with ~ 20% of the resolution 
ele ment filled by the 1700 G field an d the rest by the 500 G field. 

iBellot Rubio & Colladosl (l2003l) performed numerical simu- 
lations of Stokes profiles at 630 nm and at 1 .56 /zm, assuming a 
magnetic field filling 5% of the resolution element. For consis- 
tency with the results by Khomenko et al. (2003), they param- 
eterized the field strength distribution as a decreasing exponen- 
tial with a mean value of 250 G. After adding some noise (to 
a level of 10"^ Ic, with I^ the continuum intensity) to the syn- 
thetic profiles and inverting them, the resulting magnetic field 
strength distribution retrieved from the synthetic visible data set 
was centered at kG. On the contrary, the distribution inferred 
from the synthetic infrared spectra was close to the original one. 
They concluded that, for a sufficiently small noise level, both 
infrared and visible observations should agree and retrieve the 
initial distribution. However, and advancing some of the results 
that we will show in this paper, other effects apart from the noise 
level make the 63 nm mag netomet ry unreliable. This follows 
from the results bv lMartinez Gonzale z et al.l (|2006|) that, for the 
present observing conditions (with magnetic flux densities of ~ 
10 Mx/cm^, noise level of 6.5 x 10"^ Ic, spatial resolution from 
0.5" to 1"), the 630 nm pair of lines does not carry enough in- 
formation to retrieve the magnetic field strength. 

Apart from numerical simulations, simultaneous and co- 
spatial observations in both spectral ranges should give the final 
answer to the 1 .56 yum-630 n m controversy. The first work with 
this aim was carried out bv ISanchez Almeida et alj (l2003b us- 
ing data from the Vacuum Tower Telescope (VTT) and THEMIS 
telescopes at the Observatorio del Teide. This yielded two data 
sets in the 630 nm and 1 .56 jum lines, but with a significantly 
different spatial resolution since the THEMIS telescope was not 




Fig. 1. Can K line core image from the VTT slit jaw system. 
Bright features are good indicators of magnetic activity. The 
black square marks the area scanned for this study. The lack of 
intense bright points in the observed area indicates a very quiet 
region without intense magnetic activity. 



equipped with an image stabilization system at that time. The 
separate inversion of both data sets with a Milne-Eddington ap- 
proximation led to field strengths of about 300 G for the infrared 
lines and 1 100 G for the visible ones. More surprisingly, 25% of 
the analysed signals showed opposite polarities in the same pixel 
at both wavelength ranges. Khomenko et al. (2005) showed, us- 
ing magnetohydrodynamic simulations, that a different spatial 
resolution in both wavelengths can result in apparent opposite 
polarities at the same spatial poin t. In a subsequent publication, 
iDominguez Cerdena et al.' ('2006') degraded the spatial resolution 
of the VTT data used by Sanchez Almeida et al. (2003) to match 
that of the THEMIS observations. From a simultaneous MISMA 
inversion they also obtained an important fraction of kG fields in 
the visible and weaker fields in the infrared. Consequently, they 
concluded that the information contained in the infrared data was 
incompatible to that contained in the visible data set. They assert 
that kG and weaker magnetic fields were coexisting in the res- 
olution element, the visible being sensitive to the kG fields and 
the infrared to the weaker ones. 

Recently, the polarimeter on-board HINODE has provided 
spectro-polarimetric observations at 630 nm with the highest 
spatial resolution ever achieved with this kind of data (0.3"). 
If the occupation fraction of the magnetic elements in the reso- 
lution element would have increased at least one order of mag- 
nitude from 1" to 0.3" one would have been capable of sepa- 
rating the effects of the thermod y namics a nd the magnetic field 
strength dAsensio Ramos et aTl l2007bl) . lOrozco Suarez et akl 
(2007) perform a Milne-Eddington inversion of HINODE's data 
on a quiet Sun region. The noise level is around one order of 
magnitude higher than the one of the observations presented 
in this paper. In order to perform the inversions they assume 
one single magnetic atmosphere in the resolution element and 
some contamination of stray light. The retrieved magnetic field 
strength distribution has a peak at very weak fields, in disagree- 
ment with the previous results using the pair of spectral lines at 
630 nm. However, even if this result is compatible with the mag- 
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Fig. 2. Maps of continuum intensity and integrated absolute cir- 
cular polarization for the aligned infrared (top) and visible (bot- 
tom) data. A grid is superposed for an easier visibility of the 
co-spatiality of structures. 

netic fields obtained using the infrared lines at 1 .5 fim, the va- 
lidity of the visible spectral lines at the HINODE observational 
conditions should be tested. 

In this work we present simultaneous and co-spatial 1 .56 fim 
and 630 nm observations taken at the same telescope, under very 
similar seeing conditions and using image stabilization systems. 
Spurious effects arising from a different spatial resolution in both 
data sets, due to different observing conditions, are thus avoided. 
Prior to the analysis, noise is reduced to very low values using 
a PCA procedure to ensure that it does not bias the study of 
the compatibility of the internetwork observations at these two 
wavelengths. 

2. Observations 

On August 17, 2003, a very quiet region at disc centre was 
observed at the VTT of the Observatorio del Teide from UT 
07:34:14 to 08:26:31. A two-dimensional map was obtained by 
stepping the solar image perpendicularly to the spectrograph slit. 
At each point, the full Stokes vector was recorded at the 1 .56 yum 
and 630 nm spectral ranges simultaneously. The light coming 
from the telescope was divided using an achromatic 50% - 50% 
beam splitter, which sent half of the inco ming radiation to the 
infrared polarimeter TIP ("Collados', '1999*) and the other half to 
the visible polarimeter POLIS (Beck et al., 2005b). 

In order to study the internetwork, network regions or other 
activity areas on the solar surface were explicitly avoided with 
the aid of live Ca ii K images, which were accessible during the 
observations. The emission in the core of this spectral line is a 
good indicator of the magnetic activity (e.g., Lites et al., 1999). 
In Fig.[T]the observed region as seen in the Ca ii K line is shown. 
The black square contains the scanned area and the vertical line 
inside it represents the position of the slit at a given time. There 
were no intense bright zones in the whole scanned region, indi- 
cating that there was no significant magnetic activity. The size 
of the scanned area was 33.25" (along the slit) x 42" (along 
the scan direction). The image was stabilized using a correlation 



tracker system dBallesteros et aI.L Il996l) . which allowed accu- 
rate stepping by 0.35". The sampling of the spectral images was 
0.35" X 29.6 mA at 1.56 yum and 0.14" x 14.9 mA at 630 nm. 
The integration time at each slit position was around 27 s. The 
sHt widths corresponded to 0.35" for TIP and 0.48" for POLIS. 
The spatial resolution was estimated to be about 1.2" and 1.3" in 
the infrared and visible spectral ranges, respectively. The contin- 
uum contrast of the granulation were 1.3% (infrared) and 2.7% 
(visible). 

3. Data reduction 

Dark counts subtraction, flat-field coiTection and polarimetric 
demodulation were performed with software packages avail- 
able for both instruments. Most of the instrumental polarisa- 
tion crosstalk was re moved using cahbration optics loca ted be- 
fore the beam splitter dSchhchenmaier & CoIIadosll2002l) . For a 
complete correction, the coelostat configuration of the telescope 
needs to be modeled (C ollados , 1999; Becketal. , 2005 a). The 
residual crosstalk from Stokes I to Stokes Q, U and V was re- 
moved as: 



P' = P - yl. 



(1) 



being P or P' either Stokes Q, U or V. The symbol / accounts for 
the intensity profile. The y coefficient is obtained by forcing the 
continuum of the polarization profiles to zero since we do not 
expect continuum polarization due to the Zeeman effect. Then, 
y can be computed as: 



r = 



(2) 



being P^ and I^ the values the continuum of the Stokes Q, U or 
V profiles and of the intensity profiles, respectively. The residual 
crosstalk between Stokes Q, U and V is very difficult to remove. 
Mainly in the 630 nm data, a few percent may still remain. 

After this standard data reduction, two key items still re- 
mained to coiTect for a precise and reliable analysis: the differ- 
ential refraction of the Earth's atmosphere and the noise level. 



Differential refraction and data alignment The differential re- 
fraction of the Earth's atmosphere leads to a time-dependent spa- 
tial displacement of the solar images in different wavelengths, 
whose amount and direction depend on the combined geo- 
metrical/optical configuration of the telescope and instruments 
dReardonl l2006h . We already tried to minimise refraction effects 
during the observations by manually moving the POLIS scan 
mirror correspondingly in-between exposures. In order to correct 
for the residual differential diffraction, the continuum images of 
both spectral ranges were aligned. To that aim, the continuum in- 
tensity map from TIP data was divided into 24 small slices with 
a size of 33.25" x 1.75". The POLIS continuum image was di- 
vided into slices of 45.64" x 1.75". The correlation coefficient 
r(c/v, (iv), for integer pixel displacements, was calculated as 



r{dx, dy) 



Y^ijiAAi + d^, j + dy) - AT][Ap(i, j) - Ap] 



ylZijlArii + d,, j + dy) - Ar]2 ^i:,-,[Ap(/, ;) - Ap]^ 



,(3) 



where Ap and Ap are the average intensities of the matching 
areas between the TIP and POLIS slices, respectively, and d^ 
and df denote displacements along and perpendicular to the slit, 
respectively. The position of maximum coiTelation coefficient 
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gives the required displacement to have co-spatial granulation 
maps. These displacements inferred from the continuum images 
were then used to retrieve the co-spatial spectra at each pixel. 
The displacement perpendicular to the slit was always smaller 
than 0.6". Figure |2] shows the aligned maps. It is clear that each 
recognisable feature is located at the same position in both the 
1 .56 jjm and the 630 nm continuum intensity (or Stokes V) map. 

Noise level The noise level reached during the observations was 
2 X lO""* and 2.5 x 10""* Ic for the 1.56 jum and 630 nm data, 
respectively, with Ic being the mean continuum intensity. This 
noise level is typical for spectro-polarimetric data. The expected 
signals for internetwork fields are however so weak («: 10"^ Ic) 
that we tried to improve as much as possible the signal-to-noise 
ratio of our data. In order to reduce the noise level, we used a sta- 
tistical de-noisin g procedure based o n the Principal Component 
Analysis (PCA) (iRees & Yind.l2003h . The main idea is to sepa- 
rate the noise contribution (or residual features from the reduc- 
tion process) from the real polarimetric information as follows. 
Each one of our n observed profiles is assumed to be a 
vector of size p, the number of wavelength points. In prin- 
ciple, each observed profile would be represented by a point 
in a p-dimensional space, so that one would need p param- 
eters to fully describe it. If coiTelations exist between these 
parameters, the cloud that represents all the observed profiles 
in the p-dimensional space must be elongated in some partic- 
ular directions. If the matrix O contains the n x p observed 
profiles, these directions of maximum correlation can be ob- 
tained b};_diagonalisingJhecoiTelation matrix of the observa- 
tions JCasini & Lopez Aristel 120031) . defined as follows: 



X^O' 0, 



(4) 



where X is the p x p correlation matrix, the operator ■ indicates 
the matrix product and t denotes matrix transposition. The eigen- 
vectors of the coiTelation matrix form a suitable basis in which 
the observations can be decomposed as follows: 



C = OB, 



(5) 



being C the n x p coefficient matrix and B the p x p matrix 
containing the eigenvectors of the correlation matrix. Each C,j is 
the projection of the i-th observation on the j-th eigenvector. 

The election of this basis (the principal components) is 
not casual. The correlation matrix is constructed from solar 
profiles and the eigenv ectors may have a physical meaning 
JSkumanich & Lopez Aristel |2002|) . However, the most impor- 
tant property of the PCA decomposition for our purposes is that, 
sorting the eigenvalues in descending order by their value, their 
magnitude rapidly decreases. This means that the first few eigen- 
vectors contain the majority of the information. Under the as- 
sumption that the relevant information is contained in the eigen- 
vectors with higher eigenvalues, the PCA basis can be truncated 
and those eigenvectors with small eigenvalues can be eliminated. 
With this, one ends up with a set of p' eigenvectors that contain 
mainly useful solar information. The data can be reconstructed 
as 



Cfiit = C • B' 



(6) 



where Onu is thenx p matrix with the filtered data set, C' is the 
nx p' coefficient matrix and B' is the p x p' pseudo-basis ma- 
trix that contains only the p' first eigenvectors. The de-noising 
procedure is very easy, but the election of the new basis to repre- 
sent the data set is a tricky point. In an ideal case, one is dealing 



with spectro-polarimetric information contaminated by an addi- 
tive noise. In this situation, the PCA decomposition separates 
perfectly the noise and the real information, giving the smaller 
eigenvalues to the eigenvectors dominated by noise. In a real 
case, there exist unfortunately some other sources of noise (for 
instance, interference fringes, cross-talk from one polarisation 
state to the others, bad pixels, residual features of telluric lines on 
the polarization profiles, etc). All these contributions to the ob- 
served spectra have some kind of pattern and are not compatible 
with uncoiTelated noise. The fact that they present clear patterns 
means that sometimes those spurious signals appear mixed with 
spectro-polarimetric information with non-negligible eigenval- 
ues. Thus, the election of the new basis partly has a subjective 
character We kept only the first 16 eigenvectors for the 1.5 yum 
data and the 31 first ones for the 630 nm data. In both cases 
the eigenvalues of the rejected profiles are more than 3 orders 
of magnitude smaller than th e maximum eigenvalue. Recently, 
Asensio Ramos et alj (l2007cl) have proposed a numerical proce- 
dure to find the optimum number of eigenvectors for such kind 
of problem. Our selection of the reduced basis is in agreement 
with their calculations. 

After applying this procedure, a filtered set of spectra with 
a very low noise level is obtained. The Stokes V profiles show 
noise values of 4 x 10"^ and 7 x 10"^ Ic at 1.56 /im and 630 
nm, respectively. The noise in the linear polarization spectra is 
10""* Ic in both data sets. 

4. Analysis of polarization profiles 

4. 1 . Amplitude of circular and linear polarization 

In this section we focus on the study of directly measurable po- 
larimetric quantities: the amplitudes of circular and linear po- 
larization as well as the amplitude and area asymmetries of 
the Stokes V profiles. Interestingly, 92.6% of the observed area 
showed signals above three times the noise level, indicating that 
the whole surface is full of magnetic fields. However, in order 
to make a reliable analysis we selected those profiles that fulfill, 
simultaneously in both spectral ranges, the requirement that the 
maximum total polarization degree is larger than lOcr, with cr the 
standard deviation of the noise level at continuum wavelengths. 
This translates into 4 x 10""* for the infrared and 7 x lO""* Ic for 
the visible. In order not to bias our results, the amplitude of to- 
tal polarization Ap was computed for the line with less magnetic 
sensitivity in each spectral range (1.5652 //m and 630.15 nm) as 
the maximum of ^JQ^ + U^ + V^. This gives 56.4% of the field 
of view with total polarization signal over the lOo" threshold. 

Although the profile selection was carried out with the spec- 
tral lines with the weakest signals, we computed the ampli- 
tudes of circular and linear polarization with the lines with the 
largest magnetic sensitivity to the magnetic field (1.5648 jum 
and 630.25 nm). The amplitude of circular polarization (Ay) has 
been calculated as the maximum of Stokes V and the amplitude 
of linear polarization (Ai) as the maximum of -\JQ^ + U^. Figure 
[3] shows the resulting histograms. The circular polarization dis- 
tribution is very similar in both spectral ranges, although small 
differences exist. The visible lines detect slightly less weak sig- 
nals than the infrared ones and slightly larger amplitudes. There 
is also an interesting feature common to both spectral ranges: 
the peak of the distributions is clearly above the noise level. One 
would expect that undetected weak signals (due to the noise level 
or to the sensitivity of the det ectors) should lead to maxima of 
the distributions at noise level. lWang et alj ( ll995l) studied the in- 
ternetwork magnetism at a spatial resolution of ~ 2" using Big 
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Fig. 4. Histograms of the amplitude (left) and area (right) asymmetries of the two-lobed Stokes V profiles at 1 .5648 ;um (solid line) 
and at 630.2 nm (dashed line). 



Bear deep magnetograms. They observed a magnetic flux distri- 
bution that peaked at 5 x 10'^ Mx (~ 10 Mx/cm^), far from their 
detection threshold of lO'^ Mx (~ 2 Mx/cm^). They examined 
seeing and selection effects on their results and concluded that 
they were not responsible for the peak and the consequent drop 
towards the smallest values of the distribution down to the de- 
tection limit. We have also performed numerical simulations to 
study the solar nature of this peak above the detection limit. In 
our case, we obtain that cancellations due to the poor spatial res- 
olution produce a similar drop towards the detection limit. This 
means that the observed peak well above the noise level might 
be an observational evidence of cancellations in the resolution 
element. 

The linear polarization histograms show very different dis- 
tributions in the two spectral ranges. The visible lacks many of 
the signals that are detected in the infrared. Of importance is to 
note that the Stokes Q and U signals observed in the infrared 
present amplitudes that are of the same order of magnitude as 
those found in Stokes V. The peak of the visible distribution of 
Al is close to the noise level of 10"^ Ic, and almost no signal is 
detected for values of A/, above 10""^ Ic, just the value for which 
the infrared distribution reaches its maximum. The infrared dis- 
tribution clearly presents a peak for values much higher than the 
3cr level, as akeady found for the amplitude of circular polariza- 
tion. 

In order to shed some light into the previous re- 
sults, it is interesting to recall the q uantities introduced by 
iLandi deglTnnocenti & Landolfil (1200 4) that give an idea of the 
sensitivity of a spectral line to the circular and linear polariza- 
tion. These quantities are defined as: 



Sv 



sl 



Aref 






(7) 



where sv and sl are the circular and linear polarization sensi- 
tivity index, respectively. The symbol /l„y stands for a reference 
wavelength. Consequently, these numbers have only sense when 
comparing two different spectral lines. The symbol "g stands for 



the effective Lande factor of the transition and dc is the central 
depression in terms of the continuum intensity: 



/, - I{Ao) 



(8) 



with /(/lo) being the intensity at the minimun of the line mini- 
mum. The symbol G plays the same role as the effective Lande 
factor but for linear polarization. It is defined as: 






G^t-S, 



(9) 



where 6 is a. qu antity that depends on the quanturn numb ers of the 
transition (see iLandi deglTnnocenti & Landolfil |2004|) . In both 
the 1564.8 nm and 630.25 nm spectral lines we have 6 - Q. The 
ratios between the circular and linear polarization indices of the 
infrared and visible spectral lines are: 



isv) 



1.5/im 



(■Sv)630nm 
(■SL)l.5/jm 
(«L)630nm 



L37 
4.17. 



(10) 



This means that the sensitivity to the circular polarization is sim- 
ilar in both wavelengths but the sensitivity to the linear polariza- 
tion is 4 times larger at L5 ;um than at 630 nm. This result can 
explain the observed behaviour: the presence of ubiquitous lin- 
ear polarization signals in the infrared and the lack of them in 
the visible. 

The bottom panel of Fig.[3]shows the ratio of the circular and 
linear polarization amplitudes. The histogram of the 630.25 nm 
line shows a broad maximum near ~ 3, i.e. the majority of the 
circular polarization signals are 3 times larger than the linear 
polarization signals. The L5 /^m histogram has a peak at ~ 1, 
i.e. the linear polarization signals are as large as the circular 
ones. If we compute the theoretical ratio, sy/si, of the infrared 
and the visible lines, we obtain: 



(Sv/SL)l.5^lm 
(«v/«L)630nm 



0.34. 



(11) 



Our observational data show a ratio between the circular and 
linear polarization of the infrared relative to the visible of 0.33, 
close to the theoretical prediction. As the formulae used are 
based on the weak field approximation, this indicates that the 
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Fig. 5. The four upper panels represent the four Stokes profiles of a 630 nm observation. The four lower panels are the Stokes 
profiles of the co-spatial 1.56 fim pixel. In all the plots, diamonds represent the experimental data and solid lines the best fit from 
the two inversion procedures: the black line shows the separate inversion and red one the simultaneous one. 
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4.2. Asymmetries of the circular polarization 

In order to compute the asymmetries of the Stokes V profiles, we 
reduced the sample by selecting only those profiles that have two 
lobes in both spectral ranges. The procedure to find the ampli- 
tude and position in wavelength of the lobes is as follows. First, 
we calculated the positions where the Stokes V profile reaches 
zero. Then, we computed the maximum (or minimum) of the 
profile between two zero-crossings. With this procedure, the two 
spectral lines in the same wavelength range (e.g. 630.15 nm and 
630.25 nm) can show a different number of lobes. To select only 
those lobes that are compatible in the two spectral ranges, we 
imposed additionally the condition that the zero-crossings have 
to be the same within an error of 1 .7 km/s in the infrared data 
and 2.1 km/s in the visible data (equivalent to 3 pixels). Finally, 
only those profiles with two lobes in both spectral ranges were 
selected. With this, we ended up with 24% of the observed area. 
We used t he definition of the amplitu de {5a) and area (6A) asym- 
metries of lSolanki & Stenflol (|1984|) . The amplitude asymmetry 
is defined as: 



6a = 



(Av)h - (Ay), 

(Av)b + (Av)/ 



(12) 



where (Av)h and (Av)b represent the amplitudes of the blue and 
red lobes as defined before, while the area asymmetry is defined 
as: 



6A^ 



Ab+A/ 



(13) 



where Ab and Ar are the absolute values of the area of the blue 
and red lobes, respectively. The boundaries of the integral of 
each lobe are different for each profile. We choose the initial 
point of the integration of the blue lobe as the position in the 
blue wing closest to the lobe which has an amplitude lower than 
3 X 10"^ Ic and the ending point the zero-crossing of the profile. 
This last wavelength is the lower limit of the integration of the 
red lobe and the final one is again the wavelength point in the 
red wing closest to the profile having an amplitude lower than 
3 X 10-5 Ic. 

Figure |4] shows the histograms of the amplitude and area 
asymmetry. The amplitude asymmetry distributions for the visi- 
ble and infrared data peak at 0.15. Both are similar in shape, but 
the one corresponding to the visible lines is slig htly broader. The 
infrare d distribution is very similar to that of iKhomenko et akl 
(l2003h . who also found a peak at 0.15. The histograms of the 
area asymmetries are more symmetric than those of the ampli- 
tude asymmetries and are both centered at 0. Again, the visible 
distribution is broader than the infrared one. For comparison, the 
infrared histogram shown in Khomenko et al. (2003) peaks at 
0.07 and seems to be narrower Positive amplitude asymmetries 
and zero area asymmetries mean that our Stokes profiles have 
a higher amplitude and narrow blue lobe and a small amplitude 
and broad red lobe. This kind of asymmetric profiles are typi- 
calof active regions Hke plages and network tSolanki & Stenflol 
il984h . 



magnetic fields that give rise to the observed signals are most 
probably in the weak field regime. This means that the majority 
of the field strengths may be below x 600 G. 



5. Inversion of the full Stokes vector 

Inversion procedures are very powerful diagnostic techniques 
that allow us to retrieve the physical conditions of the solar at- 
mosphere from the information contained in the Stokes vector 
However, the smaller the signal-to-noise ratio the smaller the in- 
formation contained in the spectra. This was our main reason 
for using the PCA to reduce the noise contribution. We used the 
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jum inversion (left plot) and from the 630 nm one (right plot). 



SIR code dRuiz Cobo & del Toro IniestaL 1 19921) to carry out the 
inversions of our data sets. We adopted a two-component model 
to reproduce the observed signals in each pixel: a magnetic at- 
mosphere covering some fraction of it and a field-free one filling 
the rest of the surface. This modeling only allows us to repro- 
duce regular anti-symmetric V profiles. Consequently, we only 
used two-lobed V profiles with a signal-to-noise ratio above lOcr 
in both spectral ranges. 



The free parameters of the inversion for the field-free at- 
mosphere are: the temperature height profile (with a maxi- 
mum of 5 nodes), the line-of-sight (LOS) velocity height pro- 
file (with a maximum of 3 nodes) and the microturbulent veloc- 
ity. Concerning the magnetic component, the parameters are: the 
temperature height profile (with a maximum of 5 nodes), the mi- 
croturbulent and the LOS velocity, the magnetic field strength, 
the inclination of the magnetic field vector with respect to the 
LOS and the field azimuth. The magnetic field vector and the 
LOS velocity are set constant with height. One single macrotur- 
bulent velocity was applied to the synthetic profiles. The filling 



factor of the magnetic component is also a parameter of the in- 
version. 

As a first step, we performed separate inversions of the 
1 .56 fim and 630 nm data. As will be demonstrated below, since 
the information contained in both spectral ranges is compatible, 
we also carried out an inversion of both spectral ranges simulta- 
neously. 

The assumption of magnetic field properties and LOS veloc- 
ities constant with height prevents the recovery of the asymme- 
tries of the Stokes V profiles. Thus, profiles with strong asym- 
metries cannot be well reproduced, although other properties 
like the amplitude ratio between the two lines (for each pair of 
lines) or the width of the Stokes V lobes can be correctly fitted. 
Figure|5]shows examples of observed and best-fit profiles for the 
1 .56 fim and 630 nm lines for the separate (black line) and the 
simultaneous inversion (red line). In both cases, Stokes I and V 
are correctly reproduced. The observed linear polarization sig- 
nals are extremely low and show some crosstalk contamination 
(mainly at 630 nm). In spite of this, the shapes and the ampli- 
tudes of the profiles are nicely reproduced. 
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Fig. 8. Magnetic flux density maps retrieved from the separate inversions of the 1 .56 fim and 630 nm data. Both maps have been 
saturated to the range -50 and 50 Mx/cm^. 



6. Results 

6. 1 . Magnetic field strength distributions from independent 
inversions of tlie 1.5 ^m and 630 nm data sets 

Figure|6]shows the magnetic field strength distributions retrieved 
from the separate inversion of the infrared and visible data sets 
(solid and dashed lines, respectively). The left panel represents 
the histogra m of the magnetic field strength for a direct com- 
parison with Khomenko et al.l (l2003l) . The right panel represents 
the probability density function (P DF) for a direct comparison 
with lSanchez Almeida et all (l2003h . The PDF is the probability 
of finding a particular range of magnetic field strengths between 
B and B + dB in the field of view. It has been computed from the 
inversion results as: 



PDFiB, B + AB)^ 



Y.J(Bi,Bi + ABi) 
NAB 



(14) 



where /(B,-, B, + ABj) is the value of the magnetic filling factor 
for the magnetic field strengths that are between B, and B, + ABj. 
The symbol A^ represents the total number of profiles in the field 
of view and AB denotes the bin size. Note that the PDF takes 
into account the filling factor of the magnetic elements. 

The 1.56 jum inversions reveal that the majority of mag- 
netic fields at the internetwork have field strengths well 
below kG. This has been alread y pointed out by other works 
using the same pair o f fines JLinl 119951 : iLin & Rimmelel 
Il999t iKhomenkoetali |2003). It must be noted that the 
shape of the infrared di st ributio n resembles that pre- 
sented by iKhomenko et al.l (l2003h . On the contrary, the 
field strength distribution obtained from the visible lines 
has a peak at k G fields. This result is in agreem e nt wit h 



the works of 



iDominguez Cerdena et al. 
(I2003h. 



Socas-Navarro & Sanch ez Almeida (2002), 

(I2003b._ ISdnchez Almeida et al, 

( 12004 aiid 



Lites & Soca s-Navarrd 



iDominguez Cerdena et al.l (12006 ). The approximate height 
of formation of the infrared pair of lines is deeper than for the 
visible ones. However, this difference in height (~ 200 km) 
is not big enough to explain the extremely different values of 
the field strength. At this point, the magnetic field strength 
recovered separately from the 1 .56 //m and 630 nm inversions 
are apparently incompatible. 

In order to check the reliability of the inferred magnetic field 
distributions, we study the distributions of granular and inter- 



granular regions. The spatial resolution in our data was good 
enough to show a correlation between the continuum inten- 
sity and the bulk velocity, in the sense that dark areas in our 
data statistically correspond to downflows (intergranular lanes) 
and bright areas to upflows (granules). We select as granules 
those pixels where the continuum intensity was larger than the 
mean continuum intensity and intergranules those which present 
smaller values. Figure|7]shows the magnetic field strength distri- 
bution for granular and intergranular regions. In the infrared we 
can clearly distinguish two different distributions. In the granular 
cells, the fields are intrinsically weaker and have a distribution 
with a tail with an exponential decay. In the intergranular lanes, 
the distribution is centered at values near the equipartition field 
on the photosphere (~ 500 G). This result is expected from mag- 
ne tohydrodynamical cons ider ations and was already pointed out 
bv lLiii & Rimmelel (Il999h and lKhomenko et al.l (2003). The sce- 
nario is very different in the visible data: we found the same dis- 
tribution in granules and intergranules. This strange behaviour 
made us think that the inversion of the 630 nm lines alone might 
not be reliable. 

The magnetic filling factors infeiTed from the analysis of the 
1.56 /urn data is 1 - 2% for magnetic field strengths larger than 
300 G. Smaller field strengths are in the weak field regime and 
the magnetic field strength can not be separated from the fill- 
ing factor In the visible, the filling factor is 0.5 - 1 % for field 
strengths higher than approximately 500 G. 

6.2. Compatibility of the observations 

In this section, we present arguments about the compatibil- 
ity of our data taken at 1 .56 fim and 630 nm to be sure that 
the simultaneous inversion of both profiles at different wave- 
lengths is really trustful. Some authors working mainly with 
the 630 nm lines suggested that both spectr al ranges carry 
different information about the magnetic field tSocas-N avarrol 
2003). These conclusions are based on the systematical pres- 
ence of opposite polarities of co-spatial circular polarizatio n 
profiles at both wavelengths dSocas-Navarro & Litest |2004|) . 
ISanchez Almeida et al.l (|2003|) analysed simultaneous observa- 
tions in both spectral ranges, discovering the presence of oppo- 
site polarities in 25% of their selected profiles. Another ingredi- 
ent they point out is the higher flux density inferred from visible 
wavelengths as compared to the one recovered from the infrared. 
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Fig. 9. Magnetic field strength distributions recovered from the 
inversion of the 630 nm data set with a kG (solid line) and a 
200 G (dashed line) magnetic field strength initialization. The 
inset window represents the goodness of the fit in terms of the 
standard deviation of the difference between the observed and 
fitted Stokes V profiles. 



Fig. 10. Magnetic microturbulent velocity distributions recov- 
ered from the inversion of the 630 nm data set with a kG (solid 
line) and a 200 G (dashed line) magnetic field strength initial- 
ization. The inset window represents the difference between the 
mean temperatures of the magnetic component infeiTed from the 
weak and strong initializations. 



ISanchez Almeida et al.l (1200 3) measured an unsigned flux den- 
sity at 630 nm that was almost two times the one at 1 .56 yum. 

In our simultaneous observations, opposite polarities are not 
a common occurrence, amounting only to 3 .6% of our two-lobed 
selected profiles. Figure [8] shows the spatial distribution of the 
magnetic flux density recovered from the separate analysis of 
the infrared and visible lines, defined as: 



fB cose. 



(15) 



being 6 the inclination of the magnetic field vector with respect 
to the LOS. Both infrared and visible lines measure the same 
magnetic flux density, as shown by the strong correlation be- 
tween both maps. Only those areas of more intense magnetic 
flux density are slightly larger in the visible than in the infrared. 
Not only the value of the magnetic flux density is compatible but 
also the polarities are the same in the same pixel in both spectral 
ranges. Consequently, the conclusion can be reached that both 
spectral ranges are tracing the very same magnetic features in 
the solar surface. The following question arises: if the observa- 
tions are compatible in both spectral ranges, why does the mag- 
netic field strength distribution at both spectral ranges seem to 

be incompatible? 

iMartinez Gonzalez et al.l (l2006h showed that the 630 nm 
lines do not carry binding information about the magnetic field 
strength at internetwork areas. They showed that with typi- 
cal present internetwork observations (magnetic flux density 
~ lOMx/cm^ and a noise level of 6.5 x 10"'' Ic) the magnetic 
field strength can not be reliably recovered and separated from 
the thermodynamic parameters of the atmosphere. In order to 
check the reliability of the retrieved magnetic field distribution 
from the independent 630 nm inversions, we perform two differ- 
ent inversions with two different initializations of the code: one 
with a kG field (the one showed in Fig. |6]l and another one with 
an initial sub-kG field strength (200 G). Figure |9] shows the two 
inferred magnetic field distributions from the two different inver- 
sion procedures. The one starting with a kG field results in a field 
strength distribution with a predominance of kG magnetic fields, 
while the one initialized in the sub-kG regime is radically differ- 



ent, presenting a peak around 200 G. The inset window shows 
the mismatch (in terms of the standard deviation, cr) between the 
observational Stokes V profiles and the fits for both initializa- 
tions. Although the whole Stokes vector is well reproduced, we 
have only taken Stokes V into account to avoid our results to be 
biased by the intensity profile. All the points remain in the diag- 
onal, meaning that indistinguishable profiles are obtained from 
both different inversions. Figure[TO]shows the microturbulent ve- 
locity histograms of the magnetic component retrieved from the 
inversion of the 630 nm data set with the weak and strong ini- 
tializations. As can be seen, the initialization with a weak field 
gives rise to larger values of the microturbulent velocity. This 
is used by the inversion code to broaden the profiles since the 
magnetic field is not strong enough. The inset window in the 
same figure shows the difference between the mean tempera- 
tures of the magnetic component retrieved from the weak and 
strong initializations. The presence of different temperatures at 
the (different) heights of formation of the 630. 15 and 630.25 nm 
lines are used to modify the relative amplitudes of their Stokes 
V profiles in order to compensate for the weak value of the field 
strength. Consequently, since the magnetic field strength cannot 
be disentangled from the thermodynamics, one cannot discrimi- 
nate between the two recovered magnetic field strength distribu- 
tions shown in Fig.|9] 

According to our results, the claim bv lSocas-Navarrol (|2003|) 
that there is an observational bias produced by the difference in 
wavelength of the two spectral regions does not sustain until re- 
liable results are found for visible lines. We have been able to 
obtain sub-kG or kG fields from the same data set in the visible 
lines. Recently, Ramirez Velez, Lopez Ariste & Semel (private 
comunication) have used the Mm spectral line at 553.7 nm and 
find results that are in agreement with those obtained from the 
1 .56 //m, with a distribution of magnetic field strengths in the in- 
ternetwork well below the kG regime. The observations at 1 .56 
fim and 630 nm are compatible and a simultaneous inversion 
based on a single magnetic component does not lead to contra- 
dictions and appears to be sufficient. 
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Fig. 11. Left panel: granular (solid line) and intergranular (dashed line) magnetic field strength distributions recovered from the 
simultaneous 1.56 //m and 630 nm inversions. The dotted line is the total distribution. Right panel: distribution of the magnetic 
filling factor. 
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Fig. 12. The left panel shows the spatial distribution of the magnetic field strength retrieved from the simultaneous inversions of 
the infrared and visible data sets, the right panel represents the map of the magnetic flux density. All points set to zero were not 
analyzed. 



6.3. Simultaneous inversion of infrared and visible data sets 

When inverting the four spectral lines simultaneously, more 
physical constraints are added, leading to a larger reliability of 
our analysis. Figure |6] shows the magnetic field strength distri- 
bution (histogram and PDF representations) infeiTed from the 
simultaneous inversion of the two different spectral ranges (dot- 
ted line). The sub-kG field strengths dominate the magnetic field 
strength distribution at the internetwork. Even if the information 
of the 630 nm lines has been taken into account in the simulta- 
neous inversion, the magnetic field strength resembles substan- 
tially the one obtained with the isolated analysis of the 1 .56 yum 
lines. 

Figure \TT\ shows the magnetic field strength distribution 
on granular and intergranular regions. Granules present a his- 
togram that has a decreasing exponential shape with intrinsically 
weak fields. Intergranules, however, present a maxwellian-type 
behaviour centered on equipartition fields. The conclusion is, 
therefore, the same as that obtained from the inversion of the 
infrared pair of lines alone. The right panel of Fig. [TT]shows the 



distribution of the magnetic filling factor Most of the infeiTed 
magnetic fields fill up only 2% of the resolution element. 

The spatial distribution of the magnetic field strength is 
shown in the left panel of Fig. [12] The scale of variation of the 
magnetic field strength is ~ 2 - 5", showing patches whose size 
is similar to that found for the magnetic flux density on the right 
panel of Fig. [121 These scales are larger than the spatial resolu- 
tion and are not only due to smearing by the seeing conditions. If 
we compare the two maps, kG features are concentrated in few 
locations coiTesponding to the highest magnetic flux densities. 
On the contrary, this conclusion is not reversible, since there are 
strong magnetic flux density concentrations associated with hG 
field strengths (for instance, positions [17", 26"] and [38", 14"]). 
The spatial distribution of magnetic field strengths is really close 
to the one retrieved from the infrared inversion. Not only the 
simultaneous and the separate infrared magnetic field strength 
distributions are substantially similar, but also their spatial dis- 
tributions are in very good agreement. The spatial distribution 
of the magnetic flux density is very similar to the one retrieved 
from the separate inversion. This is an interesting point since 
while 630 nm lines are not capable of giving good results for the 
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Fig. 13. Left panel: magnetic field strength versus the inclination with respect to the LOS. Right panel: histogram of the magnetic 
flux density. All magnitudes have been recovered from the simultaneous inversion of the infrared and visible data sets. 
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Fig. 14. Left panel: magnetic flux density distributions of granules (solid line) and intergranules (dashed line). Right panel: unsigned 
magnetic flux density distributions of granules (solid line) and intergranules (dashed line). The dotted line is the total unsigned 
magnetic flux distribution. All magnitudes have been recovered from the simultaneous inversion of the infrared and visible data 
sets. 



magnetic field strength, they are good indicators of the magnetic 
flux density. 

The left panel of Fig. [T3] shows the magnetic field strength 
as a function of the inclination of the magnetic field vector It 
seems that there is a tendency for intrinsically stronger fields to 
be more vertical. Intrinsically weaker ones have a wider range 
of angles. The right panel of Fig.[T3]represents the histogram of 
the magnetic flux density. Both polarities are equally probable in 
our data set. In fact, the total flux density in the field of view is 
-0.055 Mx/cm^. This very small flux imbalance has been con- 
sidered to point towards a local dynamo as the me chanism that 
generates the mag netic fields on the internetwork dWang et al.L 
119951; iLitesl 12003 . The total unsigned flux in the field of view 
is 4.10 G, something in between the figure computed from the 
separate infrared (2.94 Mx/cm^) and visible (4.85 Mx/cm^) in- 
versions. Figure [14] shows the magnetic flux density (signed and 
unsigned) in granular and intergranular regions. There are no 
clear differences between the two distributions in granules or in- 
tergranules. It seems that the magnetic flux density remains the 
same in the whole solar surface, irrespective of the plasma mo- 
tions. 



7. Conclusions 

Simultaneous and co-spatial observations of Fe i 1 .56 //m and 
Fe I 630 nm pair of lines have been analysed. A total of 93% of 
the observed field of view has polarimetric signals above three 
times the noise level. This means that at 1" spatial resolution the 
internetwork is full of magnetic fields. 

In order to study the internetwork magnetism, the informa- 
tion of the full Stokes vector has to be taken into account. At this 
point the 1 .56 //m lines present a big advantage when compared 
to the 630 nm ones. Both pairs of lines are equally sensitive to 
the circular polarization. However, the infrared lines show con- 
spicuous linear polarization signals while the visible ones show 
these signals only in few weak features. 

The Stokes V profiles are quite asymmetric at both spectral 
ranges, having similar distributions of area and amplitude asym- 
metries. Most of the observed two-lobed profiles have an ampli- 
tude asymmetry of about 15% and an area asymmetry mainly 
zero. 
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We have performed the separate inversion of the infrared and 
visible data in order to check their compatibility. Our conclu- 
sions are; 

- The magnetic flux density recovered from the independent 
inversion is the same at both spectral ranges (both spatial dis- 
tributions have a strong correlation). Small differences can 
be produced by the different formation heights of the two 
pairs of lines. 

- There are no systematic opposite polarities for co-spatial 
Stokes V profiles in both observations. 

- The magnetic field strength distributions are very differ- 
ent. The infrared observations result in equipartition or even 
weaker fields while the visible lines show fields in the kG 
regime. 

These three points are compatible since the validity of the 
magnetic field strength distributions retrieved from the 630 nm 
data has been questioned. We conclude that the information 
of our observations is compatible in both spectral ranges. This 
means that both the 1.56 fim and 630 nm lines trace the same 
magnetic structure on the solar surface. Anyway, it does not 
mean that the infrared distribution of magnetic field strengths 
is the real one. One should keep in mind that the Zeeman effect 
is only capable of detecting magnetic fields in ^~2% of the res- 
olution element. What we really know is that our observations 
lead to reliable results and that the magnetic field strengths that 
we infer are really present on the internetwork. Ramirez Velez, 
Lopez Ariste & Semel (private comunication), with observations 
of a Mn I line in the visible that is sensitive to the magnetic field 
strength, have concluded that the magnetic field strengths on the 
internetwork are well below the kG regime. This supports the 
fact that the infrared is not biased towards weak fields. 

In order to put more physical constraints, we have performed 
the simultaneous inversion of both data sets. The inferred mag- 
netic field strength distribution is dominated by sub-kG fields. 
The occupation fraction of those fields peaks at 2% but with a 
broad distribution. Two different magnetic field strength distri- 
butions are present on granules and intergranules. The granular 
one has intrinsically weaker field strengths than the intergranu- 
lar one, which is centered on the equipartition field of the pho- 
tosphere. The magnetic flux density distribution is the same on 
granules and intergranular lanes. There is no net magnetic flux 
density in the studied field of view, implying that there is a per- 
fect cancellation of positive and negative magnetic flux density. 

The unsig ned magnetic flux density hist ograms peak at 
~ 10 Mx/cml iMartinI (Il987h and IWang et al.1 a995.) obtained 
the same value using magn etograms with 2 - 3" spatial resolu- 
tion. This is agreement with lLites & Socas-Navarrd (12004). who 
show that there is no apparent increment in the magnetic flux 
density when improving the spatial resolution from 1 to 0.6". 

Since the Zeeman effect allows us to study ~ 2% of the res- 
olution element, what happens in the rest of the surface? Is it 
fie ld-free? If not, which kind of magnetic fields can we find in 
it? lAsensio Ramos et alj (l2007ah show the first experimental ev- 
idence of magnetic flux cancellation at the internetwork. This 
supports the fact that the peak of the circular and linear polar- 
ization histograms well above the noise level is produced by 
cancellations in the resolution element. This means that there 
are still undetected magnetic fields at a spatial resolution of 1" 
(mixed polarities or extremely weak fields). The study of the in- 
ternetwork magnetism has become a very interesting field that 
deserves more efforts to know if it can play an important role in 
the so lar global magnetism as claimed by iTruiiUo Bueno et al] 
(12004 . 
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